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INTRODUCTION 

One o< tite major efforts in the development of 
gas turbines for ground power generation Is their 
adaptation to the use of coal-derived liquid fuels. 

The wide range of properties of these synthetic liquid 
fuels leads to many problems. A recent sunmary of 
coal derived liquid fuel properties and their poten- 
tial problems in turbine coadiustion is given in Ref. 1. 
The present paper focuses on the problem of Increased 
estissions. The main question is how much the oxides 
of nitrogen (NOx) emissions will be increased by the 
increased amount of organic nitrogen compounds con- 
tained in these synthetic liquids. Another question 
is whether this fuel-bound nitrogen (FBN) will have 
any significant effect upon the carbon monoxide (CO) 
emissions from synthetic fuel combustion. 

The present work was undertaken to investigate 
these questions as one part of the Critical Research 
and Technology program funded by the Department of 
Energy at the Lewis Research Center. The primary pur- 
pose of this work is to analytically determine the ef- 
fect of combustion operating conditions on the conver- 
sion of FBN to NO^. The effect of FBN and of operat- 
ing conditions on CO emissions was also examined. 
Propane-air was chosen as the combustible mixture for 
the first series of computations. The combustor model 
was assumed to be a two-stage, adiabatic, perfectly- 
stirred reactor. Chemical modeling of the oxidation 
of propane and the formation of NO^ was done using a 
fifty-seven reaction mechanism. The use of two-stage 
combustion to minimize NO^ formatlra has been reported 
by several Investigators (2 to 4).^ The greatest NOx 

iNumbers in parentheses designate references at 
end of paper. 


reduction is achieved by burning with a rich primary 
mixture and then diluting with ait for quite lean 
second-stage combustion. 

A wide range of operating conditions was used for 
these computations. Primary equivalence ratio, q>p, 
was varied from 0.3 to 1.7. Dilution air was added in 
the secondary zone, as necessary, to give a secondary 
equivalence ratio, of either O.S or 0.7. Primary 
stage residence time varied from approximately 12 to 
20 msec and secondary stage residence times of approxi- 
mately 1, 2, and 3 msec were used. Three different 
amounts of FBN were used in the computations, namely 
0.5, 1.0, and 2.0 percent by weight of the initial 
fuel. All computations were performed for a constant 
pressure of 5 stm and initial mixture temperatures 
ranged from 608 to 630 K. 

Computations were performed using the NASA Lewis 
General Chemical Kinetics Program (5), which was modi- 
to perform stirred-reactor sx>del computations. The 
nonlinear algebraic equations given in Refs. 6 and 7 
for multicomponent, chemically rate- limited homoge- 
neous reaction in a perfectly-stirred reactor were 
used. Solution of these equations by the modified 
Newton-Raphsen iteration procedure given in Ref. 6 was 
programmed and incorporated as a subroutine into the 
general chemical kinetics program. 

Results are presented for NO^ and CO concentra- 
tions as a function of iPp, Vs end both primary and 
secondary residence times, as well as the amount of 
FBN. A limited comparison of the computed results 
with experimental results obtained under another task 
in the Critical Research and Technology program is 
also given. The comparison will show the good semi- 
quantitative agreement between our theoretical model 
and Che experlmentcl results. It should be mentioned 
that neither the theoretical computations nor the ex- 
perimental work were meant to model a practical tur- 
bine combustor. 

Moreover, Che effects of initial temperature and 
pressure variation were considered of secondary impor- 
tance. These conditions were not changed in either 
Che experimental or Che theoretical work. 


CHEMISTRY OF THE COMBUSTION PROCESS 

The chemistry of nitrogen oxide formation during 
combustion has been studied extensively for many 
years. Two good reviews of this topic are given in 
Refs. 8 and 9. In this paper we shall define nitro- 
gen oxides or NOx as the sum of nitric oxide (NO) and 
nitrogen dioxide (NO2) concentrations. Although the 
mechanism of NOx ^atmation in lean hydrocarbon com- 
bustion is fairly well known, our knowledge is still 
lacking for very rich mixtures (i.e., equivalence 
ratios from about 1.2 to 2). For lean combustion NOx 
formation is explained mainly by the extended 
geldovich mechanism 


N + Oj 8 NO + 0 

(39) 

0 + Nj == NO + N 

(38) 

N -i- OH ^ NO 9- H 

(40) 


and the additional reaction; 

0 + NOj » NO + Oj (31) 

To explain experimental NO^ concentrations higher 
than computed values for rich mixtures Fenimore (10) 
suggested the direct reaction of hydrocarbon frag- 
ments with molecular nitrogen: 

CH + Nj a HCN + N (41) 

Both products of this reaction can lead to additional 
NO formation. The additional reactions of HCN, CN, 
and NHj^ species are the unknowns in the path of com- 
plete understanding of NOx formation. Jachimowskl (11) 
showed that the addition of only reaction (41) to the 
fairly well known methane oxidation mechanism improved 
rich combustion NO prediction considerably. And more 
recently, Wakelyn, Jachimowskl and Wilson (12) devel- 
oped a comprehensive propane-air combustion mechanism. 
It employs not only reaction (41) but those involving 
HCN, CN, NCO, and NH species. When applied to com- 
puting NO concentrations for stirred-reactor combus- 
tion of rich propane-sir mixtures (up to q> • 1.4), 
this mechanism gave excellent agreement with experi- 
mental results obtained in their laboratory. These 
Investigators also applied the stirred-reactor model 
and their propane-air chemical model in computations 
attempting to match experimental flame-tube emission 
data of Anderson (13). The computed and experimental 
results agreed quite satisfactorily. 

In choosing the mechanism for the present com- 
putations we have, therefore, used a modification of 
the mechanism developed in Ref. 12. One of the key 
contributions of that work is the fact that the com- 
puted NO concentration is very sensitive to the rate 
constant, k4i for reaction (41). The authors used 
this fact to derive a new value for k4^ Which gave 
the best agreement between their computed and experi- 
mental results and this improved value of k^i was 
used in the present computations. The complete mecha- 
nism used in the present computations is given in 
Table I, along with rate constants and their sources. 
The first 40 reactions plus reaction 52 were suffi- 
cient for computations at q>- values from 0.5 to 1.0. 
The full mechanism was used tor all qip values great- 
er than 1.0. 


^Reactions are numbered to correspond to numbers 
In Table I. 


SIMULATION OF FUEL- BOUND NITROGEN 

In previous analytical modeling of fuel-bound 
nitrogen conversion to NOx (4 and 13) the FBN has been 
simulated by simple species such as N, NH, NH3, HCN, 

CN. This is necessary because of Che complexity of 
the organic nitrogen compounds in synthetic fuels and 
the lack of detailed kinetic data on the pyrolysis of 
these cosq)ounds. The rationale for using one of these 
"model compounds" to represent FBN is detailed in 
Ref. 14. The basic assumption is that one of these 
simple species represents the major pyrolysis product 
of the more complicated molecule that can react to 
form NOx* The work of Ref. 14 showed chat computed 
results for NO conversion using N, NH, CN, and HCN are 
seadquantitatively similar. Inasmuch as our combustor 
sx>del is quite Idealized, we chose the simplest spe- 
cies, nitrogen atom, as our model compound for the 
FBN. It should be mentioned that using N as the model 
compound gives about 25 percent less conversion to NOx 
than using the other model compounds. This is due to 
the importance of the reverse of reaction 38, 

NO + N - 0 + H2. 

RESULTS AND DISCUSSION 

Effect of Fuel-Bound Nitrogen 

The main interest in this paper is the effect of 
FBN on both NOx CO final concentration. This ef- 
fect is shown in Fig. 1 for NOx ^ 

In each figure computed final pollutant concentration, 
with no further dilution, is plotted against Initial 
equivalence ratio, q>p, from 0.5 to 1.7 for several 
different FBN contents. For these computations pri- 
mary reaction volume was constant at 4360 cm^, the 
shortest primary-zone residence time; secondary resi- 
dence time was fixed at nominally 1 msec. These con- 
ditions were chosen to correspond to those used in Che 
experimental flame tube work which will be described 
briefly in a following section. Secondary air was 
added when necessary to give a q>g value of 0.7. 
Primary mixtures leaner than q>p “ 0.7 were not di- 
luted. It is observed in Fig. 1 chat Che minimum NOx 
concentration is always obtained at a q>p value be- 
tween 1.4 and 1.5 for the range of q>p shown. At the 
q>p value for this minimum NOx concentration Che CO 
concentration (Fig. 2) always attains its maximum val- 
ue. Although NOx always Increases as FBN is increased, 
the highest percent conversion of FBN to NOx occurs 
for 0.5 percent FBN at Che rich q>p values. For ex- 
ample, at <Pp * 1.4 the FBN conversions are approxi- 
mately 25, 19, and 15 percent for initial FBN contents 
of 0.5, 1.0, and 2.0 percent, respectively. The de- 
crease in FBN conversion with Increased FBN content 
has been observed by others (14). The effect of FBN 
on CO formation Is quite different. Figure 2 shows 
that any FBN decreases the final CO concentration. 

For q>p 1.4, the point of greatest CO concentration, 
any amount of FBN decreases the CO concentration by 
about 17 percent. 

The effect of added FBN on final NOx concentra- 
tion can be seen better in Fig. 3 where the NOx con- 
centration Is plotted against percent of FBN for sev- 
eral different primary equivalence ratios. This plot 
shows that the greatest percentage Increase of NOx 
concentration with added FBN occurs at the extremely 
lean and the very rich primary equivalence ratios. 

The absolute Increase in NOx concentration for a fixed 
addition of FBN does not change greatly with primary 
equivalence ratio. For example we can compare NOx 
concentrations for 0 and for 0.5 percent FBN In the 
following table 
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’’p 

NO^ concentration, ppm 
OX Pm 0.5X FBN 

Absolute 

change, 

ppm 

Percent 

change 

0.6 

155 

500 

345 

222 

.9 

2130 

2400 

270 

13 

1.4 

110 

340 

230 

209 


TIm abaolute changea In NOx concentration are between 
230 and 34S ppo. But thla tranalatea to a diange of 
about 200 percent at the htgheat and loweet 9p val- 
uta and a change of only 13 percent at ■ 0.9, 
where NOx fomatlon la hlgbeat without FBN preaent. 
Alao In Fig. 3, the decreaae In NOx concentration 
gradient aa percent FBN rises above 0.5 Is particular- 
ly noticeable at • 1.4, the condition for ninlmum 
NOx formation in general. 

Effect of Secondary Dilution 

The effect of additional dilution was studied by 
repeating the second stage computations for a Tg 
value of 0.5 instead of 0.7. The results for final 
NOx concentration are shown In Fig. 4 and for 00 con- 
centration In Fig. 5. All other conditions are the 
same as In Figs. 1 and 2. Final pollutant concentra- 
tion Is again plotted against 7p for various amounts 
of FBN. Comparison of Figs. 4 and 1 and Figs. 5 and 2 
show that the trends are very such the same. However, 
both pollutant concentration levels have been signif- 
icantly reduced by the added dilution. For a better 
quantitative comparison of dilution effects we have 
replotted the 1 percent FBN NOx formation curves from 
Figs. 1 and 4 together in Fig. 6(a). Likewise the 
1 percent FBN CO formation curves from Figs. 2 and 5 
are replotted together in Fig. 6(b). From these plots 
It Is seen that the greatest percentage reduction due 
to additional dilution In both NOx and CO Is at a 9p 
value close to 1.4, the value for minimum NOx* NOx ^ 
reduced by about 47 percent and CO by 57 percent due 
to the additional dilution. 

The observed concentration reductions are, of 
course, the combined effect of physical dilution of 
the primary gas and the chemical reaction in the 
secondary-stage combustion. It is of interest to com- 
pare these two effects for the present model and this 
can be done with the help of Table II. In this table 
we have listed three different values of NOx ond CO 
concentration for several <Pp values and for second- 
ary dilution to both • 0.7 and qPg • 0.5. In ad- 
dition to the secondary exhaust gas concentration we 
have listed the primary exhaust gas concentration and 
also the diluted primary exhaust gas concentration be- 
fore the secondary chemical reaction. For NOx the 
table shows that only the dilution of the primary mix- 
ture determines final concentration for the lean and 
stoichiometric (Pp values of 0.8 and 1.0. For both 
secondary dilutions the NOx concentration Is Increased 
by at most 2 percent during the secondary reaction. 

For example, at q>p - 1.0 and <Pg - 0.7 NOx concentra- 
tion increases from 2242 to 2291 ppnm during the sec- 
ondary reaction. For the rich primary mixtures, how- 
ever, there Is significant chemical reaction forming 
NOx lo the secondary stage when dilution Is to 
9, ■ 0.7. At <Pp ° 1.4, for example, the secondary 
chemical reaction Increases NOx percent from 

327 to 447 ppmm. For dilution to (Pg = 0.5 there Is, 
however, negligible chemical reaction forming NOx 
when 9p Is 1.2 and 1.4. The probable reason for 
this can be seen In Fig. 7, where secondary reaction 
temperature Is plotted as a function of 9p for both 
secondary dilutions. The same 1.0 percent FBN content 
as in Figs. 6(a) and (b) was used. The reaction tem- 


perature Is reduced by mere than 300 K for nil 9^ 
values from 0.8 to 1.7 when 9s goes from 0.7 m 0.5. 
This drop In temperature could easily quench the NOx 
formation reactions which are very taaiperaturs depen- 
dent. The situation Is quite different for CO. 

Table II shows that the final CO concentration Is 
determined mainly by secondary stage chemical reaction 
for all conditions. This Is especially true for the 
rich primary mixtures In which 00 concentration Is 
very high. The reason for the quite different behav- 
ior of CO concentration Is the lack of temperature de- 
pendence of the main CO oxidation reaction 

CO + OH « OO 2 + H (20) 

A drop of 300 K In teaq>erature, as shown In Fig. 7, 
has a negligible effect on the rate constant of this 
reaction. 

Effect of Primary and Secondary Residence Times 

Primary sons residence time was Increased by 
Changing the reaction volume by factors of 1.2 and 
1.4. Thus a typical residence time of 12 msec at the 
smallest reaction volume was changed to about 17 msec 
for the largest volume. NO significant effect or prl- 
aiary residence time was predicted for either pollutnt 
Typical examples of the variation of NOx and 00 with 
primary residence time are shown in Figs. 8 and 9. 
Pollutant concentration Is plotted against primary 
equivalence ratio for mixtures containing 1.0 percent 
fuel-bound nitrogen using three different primary 
residence times. 

All secondary-stage computationr were performed 
using three different reaction volumes. They were 
chosen to give nominal secondary re .idence time values 
of 1, 2, and 3 msec. Typical plotf of NOx concentra- 
tion versus secondary residence time with 9p as a 
parameter are shown in Figs. 10 (s) and (b). For 
Fig. 10(a) the secondary stage <Pg value Is 0.7 and 
for Fig. 10(b) It Is 0.5. Mixtures contaln'ag one 
percent FBN were used for the computations. For dilu- 
tion to 9g ■ 0.7 , we see that the NOx results are 
somewhat Sensitive to residence time variation for the 
very rich mixtures ((Pp ■* 1.4 and 1.7). There Is a 
35 percent Increase In NOx concentration as residence 
time Increases from 1 to 3 msec st 9p ~ 1.4. When 
the secondary dilution Is to <Ps * 0.5, however. 

Fig. 10(b) shows that final NOx concentration Is com- 
pletely Independent of secondary residence time. 

These results are explained by the temperature effect 
on NOx forming reactions Just discussed. There can 
only be an effect of secondary residence time on NOx 
concentration where there la chemical reaction In the 
secondary zone. It has Just been shown that NOx con- 
centration Increases In the secondary zone only for 
rich mixtures diluted to <Ps * 0.7. The effect of 
secondary residence time on CO concentration is shown 
In Figs. 11(a) and (b). The plots are exactly analo- 
gous to Figs. 10(a) and (b). In contrast to NOx, Che 
CO concentration decreases ss residence time increases 
for sll values of <Pp and both secondary dilutions. 
Moreover, the rate of decrease of CO concentration 
with residence time Is fairly independent of both pri- 
mary equivalence ratio and amount of secondary dilu- 
tion. The trends shown here, for 1.0 percent FBN In 
the fuel, are similar to those for other added amounts 
of FBN. When no added FBN Is present, the trends are 
the same especially for secondary dilution to <Ps°0.S 
The only difference Is observed for <Ps “ 0.7. For 
this dilution the NOx concentration Increases more 
sharply with Increasing secondary residence time for 
9p “ 1.4 and 1.7 than shown In Fig. 10(a). 
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CompartBOn with Experimental Data 

Exparlmantal amlaaion meaaurcmcnu are being 
Bade by G. Wolfbrandt and D. Schults at the Lewis 
Reeearch Center In another part of the CRT program. 
Their apparatus la a two-atage flame tube. A highly 
turbulent hydrocarbon*alr mixture la burned In two 
atages, with dilution In the aecond stage exactly as 
described for the present Idealized stlrrcd*reactor 
model. Both propane*alr and propane*toluenc-alr mix* 
tures are being used In this program and pyridine Is 
added to simulate fuel bound nitrogen. The purpose 
of this experimental work Is the same as for the pre* 
sent theoretical study: that Is, to determine the ef- 
fect of operating conditions on FBN conversion to NOx- 
To test whether the idealized stlrred-reactor model 
predicts the qualitative experimental trends we have 
oiade plots comparing experimental and computed re- 
sults. Figure 12 shows experimental and computed NOx 
concentration vs. ?p and Fig. 13 shows CO concentra- 
tion vs. 9>p. These results are for propane-air with 
no added FBN. The experimental polnta are the first 
data from the flame-tube experiment. For both pollut- 
ants the qualitative agreement can be seen to be quite 
good. Moreover, the quantitative agreement steadily 
Improves as we go to rich primary equivalence ratios. 
Per this latter region agreement Is within *50 per- 
cent. 

Some possible sources of discrepancy are: (1) an 
Incomplete chemical model of the cooibustion: (2) ex- 
perimental heat losses; and (3) mixing nonhomogeneity 
during the experiment. 

The present chemical model should be equivalent 
to the model of Ref. 12 In matching the data of 
Anderson (13). Ue find this to be true when we com- 
pare the results of our primary zone computations 
with the experimental results of Anderson. Our com- 
puted primary zone NO^ and CO concentrations are in 
good agreement with Anderson's results for all equiva- 
lence ratios. Therefore the observed discrepancies 
between computed and experimental second-stage pollut- 
ant concentrations ate probably due to mixing non- 
homogeneity and heat losses in the second stage com- 
bustion. One Important point is tlie fact that there 
is a finite length of mixing and chemical reaction at 
the entrance to the secondary zone for the experi- 
ment. The computational model assumes instantaneous 
mixing of the dilution air with the primary exltaust 
gases before the secondary reaction begins. Apparent- 
ly the experimental effects are less Important wlien 
the primary-zone burning is rich. 

SUMMARY OF RESULTS 

In this study we have predicted analytically the 
effect of operating conditions on pollutant fonrition 
In two-stage combustion. In particular we have 
studied the effect of these conditions on fu> , bound 
nitrogen conversion to NOx and also on CO fo.mation. 
Propane-air mixtures were assuned to burn In i two- 
stage, adiabatic, perfectly-stirred reactor. Various 
amounts of fuel-bound nitrogen, up to 2 percent by 
weight as nitrogen atoms, were used. The product 
gases of the primary stage were diluted by added air 
In tlie second stage to an overall equivalence ratio 
of cither 0.7 or 0.5. The computed results were com- 
pared with a limited amount of experimental data. 

The Important results of this work may be summarized 
as follows; 

1. Minimum NOx formation occurs In general when 
the primary zone equivalence ratio is between 1.4 and 
1.5. However, for this range of values, the CO forma- 
tion Is at Its highest level. 


2. The percentage conversion of fuel bound nitro- 
gen to NOx greatest for small percentages of FBN 
and decreases as FBN content Is raised above 0.5 per- 
cent by weight. 

3. The final concentrations of NOx and CO are 
significantly reduced by additional secondary dilution 
from 4>g - 0.7 to q>g - O.S. The extra dilution 
lowers the secondary reaction temperature by about 
300 K and thus suppresses all secondary-stage NOx re- 
actions. Therefore, the change in NOx concentration 
due to the added dilution Is strictly an effect of 
more dilution to the primary mixture. For CO, however 
the final CO concentration Is always a combined result 
of diluting the primary mixture and secondary-zone 
chemical reaction. This Is because temperature has 
very little effect on CO destruction chemistry in the 
present model. 

4. For secondary dilution to 4>g • 0.5, NOx con- 
centration is completely Independent of secondary real 
dence time. This Is because of the suppression of NOx 
reactions due to the low secondary zone reaction tem- 
peratures. For secondary dilution of rich primary 
mixtures to Tg ■ 0.7, NOx ^°Fma In the secondary zone 
due to chemical reaction. Therefore NOx concentration 
increases with increasing secondary residence time for 
these riPh primary mixtures only. Only In this region 
Is chemical reaction ImporC-enC in the secondary zone. 
The final CO concentration alwaye decreases with In- 
creasing secondary residence time because, in the pre- 
sent model, chemical reaction to destroy CO is not 
strongly affected by temperature. 

5. Final NOx ^ concentrations show negligi- 
ble dependence on primary zone residence time. 

6. Comparison of computed results with limited 
experimental data indicates that there is good agree- 
ment in trends between the Idealized theoretical model 
and the experiment. The quantitative agreement steadl 
ly improves as the primary equivalence ratio increases 
For the values of greatest Interest (4>p > 1.4) the 
agreement between theory and experiment Is within ±50 
percent. Apparently, effects of experimental heat 
transfer losses and mixing nonhomogeneity are less for 
these rich primary mixtures. These last two factors 
were not considered In the theoretical combustion 
mode 1 . 

CONCLUDING REMARKS 

This tlicoretical study reenforces the well-known 
Idea that emissions of NOx and CO can be reduced by 
two-stage combustion, burning rich in the primary ;■ ..ne 
and lean l.t the secondary. Our computations Indicate 
that modest additional secondary dilution, from equiv- 
alence ratio 0.7 to 0.5, has a strong effect in reduc- 
ing the absolute level of NOx and CO concentrations. 
Also It was found that there Is a large increase in 
NOx formation due to a small amount (0.5 percent) of 
added FBN In the fuel. Our computed emission concen- 
trations do not have a direct correlation with actual 
combustor emissions. However, the trends do indicate 
that removing only a moderate amount of the FBN from 
a real fuel during refining may not be sufficient to 
meet government emission standards. 

Future synthetic fuels contain higher amounts of 
olefins and aromatic hjrdrocarbons than current petro- 
leum fuels. The lower hydrogen/ carbon ratio of the 
synthetic fuels will also affect pollutant formation. 
Therefore, we hope to perform additional computations 
using as the fuel a mixture of toluene (C 7 Hg) and pro- 
pane, The results for NOx "bd CO emissions can then 
be compared with experimental flame tube results using 
the same fuel mixture. 
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